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We propose a simple model of supersymmetric dark matter that can explain recent results from 
PAMELA and ATIC experiments. It is based on a U{1)b-l extension of the minimal supersym- 
metric standard model. The dark matter particle is a linear combination of the U{1)b-l gaugino 
and Higgsino partners of Higgs fields that break the B — L around one TeV. The dominant mode 
of dark matter annihilation is to the lightest of the new Higgs fields, which has a mass in the GeV 
range, and its subsequent decay mainly produces taus or muons by the virtue of B — L charges. This 
light Higgs also results in Sommerfeld enhancement of the dark matter annihilation cross section, 
which can be ^ 10^. For a dark matter mass in the 1 — 2 TeV range, the model provides a good fit 
to the PAMELA data and a reasonable fit to the ATIC data. We also briefiy discuss the prospects 
of this model for direct detection experiments and the LHC. 
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One of the major problems at the interface of particle 
physics and cosmology is dark matter and its identity. 
There are various hues of evidence supporting the exis- 
tence of dark matter in the universe, and it is well estab- 
lished that particle physics can explain dark matter in 
the form of weakly interacting massive particles [l| . The 
standard scenario is thermal dark matter in which the 
dark matter relic abundance, as precisely measured by 
cosmic microwave background (CMB) experiments 
is determined from the freeze out of the dark matter 
annihilation in the early universe. Supersymmetry is a 
front runner candidate for physics beyond the standard 
model (SM). It addresses the hierarchy problem of the 
SM and has a natural dark matter candidate in the form 
of the lightest supersymmetric particle (LSP). It has been 
known that in supersymmetric models, a neutralino LSP 
can have the thermal relic abundance required for dark 
matter 

There are major experimental efforts for direct and in- 
direct detection of dark matter particle beside the grav- 
itational effect that it has on the universe. Indirect de- 
tection investigates final states (photons, anti-particles, 
neutrinos) from the annihilation of dark matter through 
astrophysical observations, while the direct detection 
probes the scattering of the dark matter particle off nu- 
clei in the dark matter detectors. The PAMELA satellite 
is an indirect experiment that has recently published re- 
sults on cosmic ray flux measurements. The data show 
an excess of positrons at energies above 10 GeV ^], while 
no excess of anti-proton flux is observed The publi- 
cation shows results up to 100 GeV and the experiment 
is expected to get data up to ~ 270 GeV. There is also 
new data from the ATIC (a balloon experiment) where 
one observes excess in e'^ + e~ spectrum with a peak 
around 600 GeV @. The backgrounds are nominal for 
both of these effects. Another balloon experiment, the 
PPB-BETS 0], also reports excess in the e"*" +e~ energy 
spectrum between 500-800 GeV. However, the excess is 



based on a few data-points that are not quite consistent 
with ATIC data. While there could be astrophysical ex- 
planations for these anomalies @, it is natural to ask 
whether they can be attributed to the effect of dark mat- 
ter annihilation in the halo. 

Model-independent analysis shows that the annihila- 
tion cross section required to explain the positron ex- 
cess exceeds the canonical value required by relic density 
'--^ 3 X 10~^^ cm^/s by at least an order of magnitude 
In the usual minimal Supergravity (mSUGRA) model, 
the situation is further complicated since dark matter 
annihilation to fermions in that model is P-wave sup- 
pressed, implying a much smaller annihilation cross sec- 
tion today as compared to that at the freeze out time. 
Even in the best case scenarios an astrophysical boost 
factor of 10'^ — 10^ is then needed to explain the positron 
excess [l^, which might be difficult to obtain based on 
the recent analysis on substructures(see e.g. [llj). More- 
over, in order to explain both the positron and anti- 
proton data, dark matter annihilation must dominantly 
produce leptons as direct product s fl^. [Tsl . There have 
been proposals [H, [l3| (also see jlSl. for new dark 
matter models in which the dark matter candidate be- 
longs to a hidden sector. Acceptable thermal relic density 
is obtained via new gauge interactions, dark matter anni- 
hilation today is enhanced via Sommerfeld effect |17i] due 
to existence of light bosons, and the annihilation mainly 
produces lepton final states. This set up can explain 
PAMELA data for dark matter mass of a few hundred 
GeV without needing a large boost factor, and ATIC data 
for larger masses |18| . 

In this paper we provide a concrete model to explain 
the recently measured anomalies in the cosmic ray. Our 
model is a simple extension of the minimal supersym- 
metric standard model (MSSM) that includes a gauged 
U{1)b-l [I^- The B ~ L extension is very well moti- 
vated since it automatically implies the existence of three 
right-handed (RH) neutrinos through which one can ex- 
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TABLE I: The B - L charges of the fields. Here Q, L and 
represent quarks, leptons, and RH neutrinos respectively; 
while H[ and H2 are the two new Higgs fields. The MSSM 
Higgs fields have zero B ~ L charges and are not shown in the 
table. 



plain the neutrino masses and mixings. It has also been 
shown that this model can explain inflation [20| . The 
model contains a new gauge boson Z', two new Higgs 
fields H[ and H2, and their supersymmetric partners. 

The B — L charge assignments are shown in Table 1. 
The superpotential is 
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where Wn is the superpotential containing RH neutri- 
nos, and fi' is the new Higgs mixing parameter. The new 
Higgs fields do not have renormalizable coupling to lep- 
ton and quark fermions as a result of their B — L charges. 
The U{1)b-l symmetry is broken by the VEV of the 
new Higgs bosons, v[ = {H[) and v'2 = {H!^)- This gives 
a mass mz' to Z' , where r7i|, = (27/4)(7^_^(wf -f- vf ), 
with gB-L being the B~L gauge coupling. We have three 
physical Higgs states, the lightest of which (p has a mass 
which is related to mz' through < rri^, cos^ 2/3', 
where tan (3' = v^jv'^ ■ For tan /?' « 1 we have <^ mz' ■ 
The other two Higgs states, $ and A, are heavy and have 
masses comparable to mz' ■ Note that unlike the MSSM 
case, radiative corrections to the Higgs quartic coupling 
do not lift because H'l and H2 are not coupled to 
fermions. Note that although the particle content in our 
model is similar to the one in [2l[, our set up is different. 
Assuming that supersymmetric particles in the MSSM 
sector are heavier than those in the U(1)b-l sector, the 
dark matter in this model arises from the new sector. 
The dark matter particle, denoted by Xi, is the lightest 
of the three new neutralinos Xi (not to be confused with 
the MSSM neutralinos which we do not discuss in this pa- 
per). It is a linear combination of the U{1)b-l gaugino 
Z' and the two Higgsinos H( , • The dark matter ther- 
mal relic abundance is dictated by the annihilation of the 
lightest neutralino Xi into a pair of (f) via s-channel ex- 
change of (f> and <i>, and t-channel exchange of Xi- (There 
is also subdominant annihilation processes to // via s- 
channel Z' exchange, and t-channel sfermion exchange.) 
The annihilation to 4'4' is not P-wave suppressed since 
the final state particles are bosons. Hence the (perturba- 
tive) annihilation cross-section does not change between 
the time of freeze-out and the present time. 

The (j) subsequently decays into fermion-antifermion 
pairs via a one- loop diagram containing two Z' bosons. 
The decay rate is given by: 
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where and Qf are the B — L charges of (p and the 
final state fermion respectively, mf is the fermion mass, 
and Cf denotes color factor. (The Higgs can also decay 
into four fermion final states via two virtual Z' , but this 
decay mode is suppressed by two orders of magnitude 
compared to the two fermion final states.) Decays to 
neutrinos are much suppressed because of the negligible 
neutrino masses. Since B — L charge of leptons is three 
times larger than that of quarks, the leptonic branching 
ratio is about 27 times larger than that for quarks of 
comparable mass. We note that can be controlled 
by the VEVs of the new Higgs fields and for comparable 
VEVs, i.e. for tan/3' « 1, it can be very small without 
any tuning on the soft masses in the Higgs sector. We can 
choose this mass to be between 0(1) GeV and 10 GeV. 
For 2'mr < < 2mh the dominant decay mode is to 
T^r+ final state. If is slightly less than 2mr, 4> can 
decay either to cc or with comparable branching 

ratios. It is possible to reduce the (p mass further to 
be below 2mc, and make n~ final state the dominant 
decay mode. 

The annihilation cross section at the present time 
has Sommerfeld enhancement as a result of the attrac- 
tive force between dark matter particles due to the 
light Higgs boson exchange. The Higgs coupling to 
dark matter, h(f>XiXii leads to an attractive potential 
V{r) = — a(e^™'''''/r) in the non-relativistic limit, where 
a = h'^/ir. Note that a is larger than the usual defini- 
tion of fine structure constant because of the Majorana 
nature of Xi- Since the neutralinos are traveling with 
non-relativistic speed v 10~'^c today, the Sommerfeld 
effect is much more important now than at the time of 
freeze out. The Sommerfeld enhancement in our model 
can be ^ lO'^. We can therefore explain the PAMELA 
data without requiring any boost factor. 

To show that our model can explain the PAMELA 
data, we pick random parameters and generate point 
models. We then calculate the relic density and the Som- 
merfeld enhancement factor for each of these models. We 
use reasonable values for the parameters, i.e. tan/?' « 1, 
mz' > 1.5 TeV, ^' 0.5 - 1.5 TeV, soft masses for the 
Higgs fields ttih' — 200 — 600 GeV, and soft gaugino 
mass A%, = 20o'- 600 GeV. We use qb-l 0.45, which 
is in concordance with unification of the gauge couplings 
(we need to use a normalization factor \/3/2 for unifi- 
cation). The Z' mass used in the calculation obeys the 
LEP and the Tevatron bounds [l^, for our charge 
assignments. 

In Figure [l] we show the relic density and the dark 
matter mass for different points. The horizontal band 
shows the acceptable range for relic density according to 
the latest CMB data 0]. In Figure [2] we show the possi- 
ble enhancement that can be obtained for these points in 
term of = / aM^o . Note that for points that satisfy 
the relic density constraint thereare many that have en- 
hancement factor > 10'^, corresponding to = 0.153 to 
0.157, within the whole range of A/^o shown in Figure[TJ 
In addition, we note that the lifetime of (j> for these points 
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FIG. 1: We show the rehc density and the neutrahno mass for 
model points generated by varying the parameters mentioned 
in the text. 




FIG. 2: We show relic density as a function of e^. We show 
different ranges for the Sommerfeld enhancement factor R by 
shades. 



is found to be ~ 10"'* — 10~^ seconds from Eq. (2). 
This is short enough to escape the tightest bounds from 
big bang nucleosynthesis (BBN) 

We select models that satisfy the dark matter relic 
density. We then use DarkSUSY-5.0.2 [251 to calculate 
the positron flux from dark matter annihilation. Note 
that each pair annihilation in our model produces 2 0's 
that yield four fermions upon their decay. For this rea- 
son, we generally need a heavier neutralino compared to 
models in which the pair annihilation directly produces 
two fermions. We normalize the positron fraction by a 
factor fcf, = 1.11 [26]. There are theoretical uncertain- 
ties in the positron cosmic ray flux calculation due to the 
assumptions about the dark matter halo profile and the 
cosmic ray propagation model. Here we assume a NFW 
profile [23| and MED parameters for the propagation as 
defined in [28|. 

In Figure [31 we show our fit to the PAMELA data for 
M^o = 1,1.5 and 2 TeV, for t"t+ final state case. We 
see that the fit is very good for a neutralino mass around 
1.5 TeV. We have chosen a point where the enhancement 
factor is 10^. In general, for larger enhancement fac- 
tors we can fit the data with a larger neutralino mass. 
For an enhancement factor of 10"*, we can obtain a good 
fit for a dark matter mass around 5 TeV. On the other 
hand, smaller masses require a smaller enhancement fac- 
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FIG. 3: We show a fit to the PAMELA data when the 
decays mostly to taus for neutralino masses to be 1, 1.5, and 
2 TeV (from top to bottom), with enhancement factor 10^. 
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FIG. 4: Same as Fig. |3]but for the case of final state muons. 

tor. However, for this tau case the fit with a smaller 
mass is not good because the spectrum is too soft. The 
antiproton data is still satisfied with this large enhance- 
ment factor since the leptonic branching ratio is much 
larger (~ 27) compared to the quark branching ratios. 

In Figure|4]we show the plot with ^"/i"*" final states. In 
order for the Higgs to mainly decay to muons, we need to 
make smaller than twice the charm mass. However, 
the lifetime of the Higgs becomes larger than a second in 
this case, which will be problematic for BBN [24]. The 
Higgs lifetime can be reduced if we increase the B ~ L 
charges by a factor of two. On the other hand, the 
positron spectrum in the muon case is harder than that 
in the tau case, and we can also fit the PAMELA data 
with a 500 GeV neutralino mass, and an enhancement 
factor of around 100. 

In Figure [5] we show the fit to the ATIC data by using 
muon and tau final states. We normalize the background 
to fit with the data at smaller energies. We again use 
enhancement factor of lO"^. We see that the muon final 
states give a better fit than the tau final states. However, 
we can see that the simultaneous fit to both ATIC and 
PAMELA is not satisfactory. In particular, it is more 
difficult to fit the ATIC data with our model. If we use a 
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FIG. 5: We show a fit to the ATIC data when </!> decays to 
muons mostly (sohd hne) and taus mostly (dashed line) for 
the lightest neutralino masses to be 1, 1.5 and 2 TeV 

large enhancement factor, the spectrum at smaller ener- 
gies will also be lifted up. From the ATIC data plot [1], 
there seems to be some missing signals at some energies 
which results in the jagged curve. This can also indi- 
cate more than one source for the excess. In any case, 
we need to wait for verification from future experiments. 
It is interesting to point out that there is a proposal for 
an upgraded version of ATIC, called ECAL [23], which 
should have much improved background rejection power 
and higher resolution. We will be able identify the model 
parameter space by combining results from all these ex- 
periments. 

We would also like to comment on other phenomeno- 



logical aspects of the model. The leading order interac- 
tion of dark matter particle Xi with quarks is via squark 
exchange in the i-channel. The new Higgs fields couple 
to the quarks at one-loop level, and hence interactions 
via Higgs exchange are suppressed (even after Sommer- 
fled enhancement is taken into account). As a result, 
the cross section for spin-independent interactions is also 
very small, well below 10"^" pb, and hence beyond the 
reach of direct detection experiments. We note that since 
left and right quarks have the same B ~ L charge, and 
x5 is a Majorana particle, there will be strictly no spin- 
dependent interactions between dark matter and ordi- 
nary matter in this model. The model however has a 
great potential to be observed with the Fermi Satellite 
experiment. The Sommerfeld enhancement would still 
be responsible for giving rise to a higher rate of photons 
in the cosmic gamma ray background. 

At the LHC, the Z' can be produced. However, the 
new light Higgs (/> will decay outside of the detector, hence 
will be missed, because of its long life time 10~'' sec). 
Thus we have another source of missing energy signal in 
this model. We also note that there are 7 neutralinos 
in this model, compared to four in the MSSM, while the 
number of charginos is still two. Therefore, using the end 
point analysis [30|, one can find many neutral states. 
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